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ABSTRACT 


The  problem  of  determining  whether  a  line  of  bearing 
measured  by  a  local  surface-based  sensor  coincides  with  a 
threat  whose  position  has  been  previously  estimated  by  an 
ocean  surveillance  sensor  is  addressed.  Uncertainties  in 
the  position  estimate  of  the  threat/  in  the  bearing  error 
and  in  the  position  estimate  of  the  sensor  are  considered 
in  measuring  the  probability  that  the  threat  lies  on  a  given 
bearing  from  the  sensor.  A  TI-59  calculator  program  is 
developed  which  calculates  this  likelihood  when  the  threat 
location  density  can  be  assumed  to  be  bivariate  normal. 
Computations  required  when  significant  time  has  elapsed 
since  the  original  estimate  of  threat  location  when  the  den¬ 
sity  can  no  longer  be  considered  bivariate  normal  are 
discussed. 
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I .  INTRODUCTION 


As  long  range  surface  launched  weapons  systems  continue 
to  be  introduced  into  the  fleet,  the  operational  commander 
is  increasingly  faced  with  the  problem  of  being  able  to 
launch  weapons  at  targets  located  beyond  the  horizon. 
Successful  employment  of  such  weapons  is  dependant  not  only 
on  the  ability  to  detect,  classify  and  localize  targets  at 
considerable  distances,  but  also  on  the  ability  to  distinguish 
the  true  target  from  a  potentially  larger  field  of  false 
targets.  While  long  range  ocean  surveillance  sensors  may  be 
of  assistance  in  the  identification  and  localization  of  tar- 
crets,  the  information  provided  may  not  be  refined  sufficiently 
to  permit  effective  targeting  of  long  range  weapons  on  that 
basis  alone.  The  on-scene  commander  must  in  general  rely  on 
additional  data  on  target  location  gathered  locally  and  close 
to  the  time  of  weapons  launch  for  accurate  targeting.  Thus, 
he  must  still  be  able  to  detect  and  track  the  desired  target 
and  be  able  to  distinguish  it  from  other  targets  within 
range  of  his  sensors. 

The  procedures  developed  in  this  paper  are  designed  to  be 
of  assistance  in  addressing  the  last  of  these  problems.  They 
are  applicable  when  the  information  available  is  an  error 
ellipse  around  a  threat  location  estimated  by  an  ocean  sur¬ 
veillance  sensor  and  bearings  only  data  generated  by  a  local 
surf ace-based  sensor.  The  question  of  whether  a  target 


detected  by  the  local  sensor  is  the  same  as  that  whose 
estimated  position  was  provided  by  an  external  sensor  can 
only  be  addressed  if  information  is  available  on  the  loca¬ 
tions  and  tracks  of  all  possible  targets  within  range  of  the 
local  sensor.  Since  such  data  is  generally  not  available, 
this  paper  does  not  attempt  to  answer  that  question,  but 
rather  develops  a  method  by  which  bearing  information  from 
different  sensors  can  be  compared  as  to  the  likelihood  of 
each  bearing  being  associated  with  the  threat  identified 
previously.  It  is  envisioned  that  these  likelihoods  can 
be  then  used  to  induce  an  ordering  among  bearing  data  gathered 
by  different  sensors  or,  conceivably,  conflicting  data 
gathered  by  one  sensor.  The  ordering  would  be  based  on  the 
likelihood  that  each  bearing  will  contribute  to  refining  the 
original  estimate  of  the  location  of  the  target  of  interest. 
This  information  could  be  applied  in  a  number  of  ways:  as 
a  guide  to  allocation  of  more  capable  sensor  resources  for 
purposes  of  obtaining  targeting  information;  as  a  guide  for 
allocation  of  weapons  against  more  than  one  threat;  as  a 
means  of  pre-processing  data  before  entering  it  into  a  target 
motion  model,  thereby  reducing  the  chance  of  introducing 
unrelated  data. 

To  determine  the  likelihood  that  a  given  line  of  bearing 
and  the  threat  coincide,  consideration  was  given  to  the 
uncertainties  inherent  in  estimation  of  target  position, 
in  the  measurement  of  bearings  by  a  particular  sensor  and 
in  estimation  of  sensor  location.  It  is  assumed  that  at  some 
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time  t  ,  an  ocean  surveillance  sensor  detects  a  threat  whose 
o 

position  is  estimated  to  be  within  an  elliptical  region  with 
*  100%  certainty.  The  estimated  position  data  sure  received 
and  converted  by  the  on-scene  commander  into  a  probability 
distribution  described  by  a  truncated  bivariate  normal  density 
function. 

It  is  further  assumed  that  the  standard  error  o.  charac- 

3 

teristic  of  the  local  sensor  is  known.  The  sensor  bearing  6 
with  bearing  error  0,  is  then  projected  from  the  sensor 
position  through  the  threat  density  function. 

Since  sensor  position  relative  to  the  target  may  itself 

be  subject  to  navigation  error,  the  uncertainty  is  introduced 

as  a  truncated  circular  bivariate  normal  distribution  centered 

at  location  (u  ,v  )  with  standard  deviation  0. 
o  o 

A  TI-59  calculator  program  is  developed  which  estimates 
the  likelihood  that  the  threat  identified  by  an  external 
sensor  lies  along  bearing  B,  given  the  threat  distribution, 
the  bearing  error,  and  sensor  position  distribution  relative 
to  the  threat. 

The  theoretical  basis  for  this  calculation  is  presented 
in  Chapter  II.  The  algorithms  used  in  designing  the  calcu¬ 
lator  program  cure  described  in  Chapter  III.  A  program  listing 
and  verbal  flow  are  provided  in  Appendix  A  along  with  instruc¬ 
tions  for  the  user.  Appendix  B  contains  a  verbal  flow  of  a 
program  designed  for  use  when  considerable  time  has  elapsed 
since  the  initial  estimate  of  the  location  of  a  moving 


target. 


II .  THEORETICAL  BASES 


The  general  approach  to  determining  the  likelihood  that 
a  measured  line  of  bearing  3  is  the  true  bearing  from  the 
sensor  to  the  threat  identified  and  localized  by  an  external 
sensor  is  discussed  in  this  chapter.  Calculations  required 
when  using  threat  position  information  both  as  initially 


generated  by  the  external  sensor  at  time  zero,  tQ,  and  as 


distorted  to  account  for  an  intervening  time  late,  t^,  are 
discussed.  Uncertainties  in  bearing  measurement  and  sensor 
position  are  included. 

Initially,  assume  that  sensor  position  is  known  with 


certainty-  Let  3^  be  the  true  bearing  of  the  threat  from  the 


sensor.  Since  threat  location  is  uncertain,  3„,  is  a  random 


quantity  with  probability  density  function  f_ 


C3t) .  Let 


=TRUE 

f  { 3;  3-,)  ds  be  the  probability  that  the  errors  in  bearing 


measurement  are  such  as  to  give  rise  to  a  bearing  on  the 


threat  in  the  interval  d3  about  the  observed  value  s  when 
the  true  bearing  is  ST.  Then  the  likelihood  of  observing  a 
bearing  3  is: 


f  *  (6)  = 


all  3, 


(s; B„)  f.  CsT) 

1  STRUE 


{1} 


The  density  f,  {3»}  is  determined  by  the  probability 
2  TRUE 

density  function  of  threat  location  at  time  t,  f„  v(x,y;t>. 


Given  the  probability  density  function  of  threat  location 
at  time  t  and  assuming  that  bearing  errors  are  normally  dis¬ 
tributed  with  mean  zero,  consideration  is  limited  to  the 
probability  that  the  observed  bearing  is  the  true  bearing 
of  the  threat,  given  the  threat  is  contained  in  a  planar 
region  E  acid  the  true  bearing  Bm  is  in  an  interval  about  3 

"  i. 

with  upper  and  lower  bounds  =  and  S_ .  The  region  E  is 
selected  to  be  the  minimum  area  planar  region  which  contains 
the  threat  with  a  specified  high  probability  pn .  The  inter¬ 
val  ( 3- * 3„)  is  selected  so  that,  for  all  the  possible  values 
of  3,p  contained  in  the  interval,  3  is  within  am  interval  of 
specified  high  probability  p2  around  3^.  The  fan  (3^,5^) 
is  symmetric  about  3  so  that  it  represents  the  minimum  area 
region  which  meets  the  amove  criterion.  The  likelihood  of 
observing  3  when  the  threat  is  in  E  amd  sT  is  in  (3^,3^)  is 
determined  by  integrating  expression  (1)  above,  over  a  region 
defined  in  the  manner  of  the  shaded  area  of  Figure  1(a). 

If  sensor  position  is  not  certain,  assume  that  it  is 

distributed  in  accordance  with  a  circular  bivariate  normal 

distribution  f „  „{u,v).  Again,  consideration  is  limited  to 

u,  v 

determining  the  likelihood  o*  observing  a  bearing  3  given 
that  the  threat  is  contained  in  a  planar  region,  E,  3^  is  in 
( ST ,  ST7) ,  amd  the  sensor  is  contained  in  a  planar  region  C. 

The  region  C  is  selected  as  the  minimum  area  pi  an  air  region 
which  contains  the  sensor  with  a  high  probability  -  The 
region  of  E  over  which  f*(3)  is  evaluated  expands  as  the 
shaded  region  of  Figure  1 (b) . 


(b)  Sensor  Po 


FIGURE  1 .  REGION 
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Since  the  measurement  errors  involved  in  estimating  the 
threat  position,  the  sensor  position  and  the  bearing  angle 
arise  from  different  measurement  procedures,  it  is  reason¬ 
able  to  assume  that  errors  are  independent.  The  probability 
densities  required  to  perform  the  above  calculations  can  be 
estimated  as  follows. 

At  time  t  an  ocean  surveillance  sensor  estimates  the 
o 

position  of  the  threat  to  be  located  within  an  elliptical 
area  characterized  by  the  parameter  set  E  =  {X,Y,9,A,B}  with 
confidence  p^  *  100%.  The  elements  of  the  parameter  set  E 
are  respectively:  X  the  latitude  of  the  estimated  threat 
position,  Y  the  longitude  of  estimated  threat  position, 
e  the  orientation  of  the  major  axis  from  true  North,  A  the 
length  of  the  semi-major  axis,  and  B  the  length  of  the  semi¬ 
minor  axis.  Since  the  measurement  errors  in  determining 
threat  position  are  generally  assumed  to  be  normally  distributed, 
the  ellipse  characterized  by  E  represents  the  minimal  area 
p-^  x  100%  confidence  region  about  the  mean  (X,Y)  .  Treating 
this  ellipse  as  the  p^  probability  region  of  a  bivariate 
normal  distribution,  a  density  function  for  the  threat  posi¬ 
tion  at  time  tQ  can  be  estimated.  For  convenience,  locate 
the  origin  of  a  rectangular  coordinate  system  at  the  center 
of  the  ellipse,  (X,Y) ,  with  positive  x-axis  located  along 
the  major-axis  of  the  ellipse  at  a  bearing  9  from  true  North. 
Assume  a  flat  earth  in  the  region  of  interest.  Let  t  =  0. 

The  mean  of  the  threat  position  density  fv  v(x,y;0)  is  then 
the  point  (0,0).  The  variances  in  the  X  and  Y  directions  can 
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be  derived  from  the  fact  that  the  region  with  minimal  area 


which  contains  the  threat  with  probability  p^  is  a  k-sigma 
ellipse  where  k  is  determined  from  the  relationship: 

.k2 

P  [threat  located  in  k--sigma  ellipse]  =  1  -  e 


Thus, 


ox2  =  (A/k) 2 

and 

<Jy2  =  (B/k)2 


where 


k  =  /-2  In (1  -  p^  . 

X  and  Y  are  assumed  to  be  independent. 

If  the  course  and  speed  of  the  threat  are  known  with 

certainty  to  be  i| >  and  s  respectively,  the  probability  density 

of  the  threat  position  at  time  late  tL  can  be  shown  to  be 

again  a  bivariate  normal  with  mean  (st^cos (8--.fi)  ,stLsin(9-^) ) 

2  2 

and  variances  ax  ?  ay  •  The  P]_  Probability  region  of  the 
density  at  time  t  would  then  be  an  ellipse  congruent  to  that 
characterized  by  the  set  E  above  but  centered  at  the  point 


In  some  cases,  the  motion  of  a  submarine  on  patrol  in  a 
large  area  can  be  characterized,  when  t^  is  large,  by  an 
expansion  of  the  probability  area  with  time  at  a  rate  D. 

The  result  in  such  a  case  is  that  fv  ^(x^.'t-)  is  still 

A/  x  u 

2 

bivariate  normal  with  mean  (0,0)  but  with  variances  ax  +  DtL 
2 

and  ov  +  Dt_  . 

x  Lt 

When  the  motion  of  the  threat  cannot  be  described  in  the 
above  manner,  and  the  course  and  speed  are  not  known,  but 
assumed  to  be  distributed  according  to  the  densities  ff(4») 
and  fg(s),  determining  the  probability  density  fx  ^Cxjyjt^) 
is  a  considerably  more  complex  problem.  Let  (xQ  (s,  y)  ,y  (s,  *41 ) ) 
be  the  coordinates  of  the  point  at  which  the  threat  would 
have  to  be  located  at  time  zero  in  order  to  reach  the  point 
( Xl , y l )  at  time  tL  if  the  threat  speed  were  s  and  course 
Thus,  xQ{s,uf)  =  xL~stLcos  and  yQ[s,-p)  =  yL~stLsin (0-iJ>) 

(Figure  3) .  Then, 

®  360 

fx^Y(x,y;tL)  =  /  /  f X , y  (x0  ( s ,  $ )  , yo  ( s ,  ’J' )  r  tQ ] 

•  ff(ii/)fg(s)  d’^ids  [Ref.  1]  . 

This  density  is  no  longer  normal.  In  the  special  case  when 
2  2 

A  ~  B,  i.e.,  ax  -  Oy  ,  'I/  has  a  uniform  distribution  over 
the  interval  (0°,360°)  and  S  is  known  with  certainty,  as 
shown  in  Reference  2,  the  density  changes  with  time  as  in 
Figure  4. 
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TIME  LATE  THREAT  DENSITY,  SPEED  KNOWN 
COURSE  UNIFORM  (0°,360°) 


The  distribution  of  the  bearing  error  measured  by  the 

local  sensor  can  be  estimated  if  the  standard  error  of  the 

sensor,  a is  known.  The  bearing  error  is  then  assumed  to 

2 

be  normally  distributed  with  mean  zero  and  variance  . 

D 

Assume  that  the  estimated  sensor  location  (u  ,vj  is 

o  o 

accurate  to  within  R  nautical  miles  with  p^  x  100%  confidence. 

The  density  fy  V(u,v)  of  sensor  location  can  be  assumed  to  be 

a  circular  bivariate  normal  with  mean  (u  ,v  )  and  variance 

o  o 

2  2  2 

o  =  R  /k  .  The  value  of  k  is  determined  from  the  relation¬ 
ship: 


P((u,v)  contained  in  k-sigma  region] 


where  the  probability  on  the  left  is  p^  in  this  cas  . 

Since  the  evaluation  of  f  * ( 3 )  considering  sensor  position 
density  and  bearing  error  density  does  not  generalize  to  a 
closed  form,  algorithms  are  developed  in  the  remainder  of  this 
paper  for  estimation  of  the  probability  (likelihood)  that  the  threat 
identified  by  an  external  sensor  lies  on  a  line  of  bearing 
measured  by  a  local  sensor  given  that  the  threat  is  in  region 
E,  the  sensor  is  in  region  C  and  the  true  line  of  bearing 
lies  within  (3^,8^. 
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III.  ALGORITHMS 


As  indicated  in  the  previous  chapter#  the  variety  of 
geometrical  situations  which  can  arise  depending  on  the 
location  of  the  sensor  relative  to  the  estimated  threat  loca¬ 
tion  precluded  development  of  a  generalized  analytical  pro¬ 
cedure.  Rather,  algorithms  are  developed  in  this  paper  for 
numerically  evaluating  f'(3).  The  procedures  applicable  to 
a  bivariate  normal  threat  location  density  have  been  imple¬ 
mented  on  a  TI-59  calculator.  Appendix  A  contains  a  listing 
of  that  program.  The  calculations  required  when  the  time 
late  threat  location  density  is  no  longer  normal  exceeded 
the  available  program  capacity  of  the  TI-59  and  therefore 
have  not  been  implemented.  A  detailed  verbal  flow  is  provided 
at  Appendix  B  for  future  implementation  on  a  larger  machine. 
The  algorithms  used  in  both  situations  are  described  in  this 
chapter. 

A.  BIVARIATE  NORMAL  THREAT  LOCATION  DENSITY 

This  case  includes  situations  (1)  where  the  time  elapsed 
since  generation  of  the  initial  error  ellipse  by  the  ocean 
surveillance  sensor  is  negligible,  (2)  where  the  motion  of 
the  threat  can  be  assumed  to  be  random  in  the  manner  described 
above,  and  (3)  where  course  and  speed  of  the  threat  are 
assumed  to  be  known  with  certainty.  With  appropriate  modi¬ 
fications  to  the  input  data,  all  three  of  these  situations 
can  be  addressed  using  the  program  contained  in  Appendix  A. 
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In  situation  (1}  the  data  entered  are  the  parameters  of  the 
ellipse  as  generated  by  the  ocean  surveillance  sensor.  In 
(2)  ,  the  location  and  orientation  from  North  of  the  ellipse 
is  the  same  as  originally  generated,  but  the  size  of  the 
ellipse  expands  at  some  constant  rate  of  area  per  unit  time 
which  must  be  estimated  by  the  user.  This  rate  D,  times  the 
elapsed  time,  t^,  yields  the  factor  which  must  be  added  to 
the  semi-major  and  semi-minor  axes  of  the  original  ellipse. 
That  is,  if  the  original  error  ellipse  is  a  p1  *  100%  confi¬ 
dence  ellipse,  the  semi-major  and  semi -minor  axes  of  the 
diffused  ellipse  will  be  input  as  A'  and  B*  respectively: 


In  situation  (3)  the  dimensions  and  orientation  of  the  time 
late  ellipse  are  the  same  as  those  of  the  original  error 
ellipse,  but  the  center  of  the  ellipse  is  displaced  from 
its  original  position  by  the  known  velocity  times  elapsed 
time.  The  updated  position  of  the  error  ellipse  is  treated 
as  the  origin  of  the  rectangular  coordinate  system  for  this 
situation  and  all  linear  measurements  are  made  relative  to 
this  system.  All  angular  measurements  are  made  from  true 
North. 

Estimation  of  the  likelihood  that  the  threat  lies  on  a 
bearing  from  the  local  sensor  given  that  the  threat  is  located 


in  a  100%  confidence  ellipse  and  the  true  bearing  lies 
within  the  bounds  (3^,3^)  with  confidence  x  100%  proceeds 
as  follows: 

1.  Estimate  parameters  of  the  bivariate  normal  density 

2  2 


lX,Y 


(x,y)  of  threat  location: 


X 


SY  "  °'  °X  ' 


2.  Determine  3_  and  3„  such  that  an  interval  of  length 

1*  u 

2k, a.  centered  on  either  3T  or  3„  would  contain  s,  the 

P  D  Li  U 

measured  bearing,  with  probability  d_:  3t  =  3  -  k_o„, 

*  6  XJ  5  P 

8  it  =  3  +  kS°S* 

3.  Determine  sensor  location  coordinates  relative  to 
the  origin  of  the  threat  ellipse. 

4.  Subdivide  the  angular  interval  (3L,3y)  into  2n  sub¬ 
intervals.  Each  subinterval  k  intersects  the  ellipse  in  a 
strip  with  average  width  which  corresponds  to  ag^. 

5.  At  the  midpoint  of  each  subinterval  k  determine  the 
equation  of  the  line  through  the  sensor  position  at  the  true 
bearing  from  North. 

Y-c 

6.  Let  the  equation  of  the  line  of  bearing  3^  be  X  =  -jp . 
Then  the  plane  perpendicular  to  the  xy-plane  which  contains 
this  line  intersects  the  bivariate  normal  threat  density  in 

a  curve  whose  equation  is 


-l(.(y-rg  )±  +zL) 
2k  2  2  2} 
ni  vj o »» 


g(y)  = 


2*axcx 


(2) 


V-Q 

found  by  making  the  substitution  X  =  -p-  in  the  density 
function  f^  y(x,y).  It  will  prove  convenient  to  expand  the 
right  side  of  (2)  as  follows  [Ref.  3]: 
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For  computational  purposes,  assume  that  the  width 
of  the  region  of  the  ellipse  cut  out  by  the  angular  subinter¬ 
val  around  3^  will  be  nearly  constant  through  the  ellipse. 

Let  the  points  of  intersection  of  the  line  of  bearing 
with  the  threat  ellipse  be  (Xj^Yj^)  and  ^Xk2/~k2^  *  Tiien 
approximate  the  volume  of  the  normal  density  over  this  region 
by  the  absolute  value  of  the  product  of  the  area  under  the 
curve  g  between  and  Yk2  and  (Figure  5} .  Observe  that 
the  term  in  brackets  in  equation  (3)  above  is  the  density 
function  of  a  univariate  normal  random  variable  with  mean 


and  variance 


2  2  2 

P  °X  °Y 

2  2  2  * 
Cy  +  m  ox 


The  term  preceding  the  brackets  in  (3)  is  the  slope  of  the 

line  of  bearing,  m,  times  the  density  of  a  univariate  normal 

2  2  2 

random  variable  with  mean  zero  and  variance  +m  ox 


evaluated  at  c 


Thus,  evaluating  the  expression 


Y 


k2 


g(y)  ay  I 


'kl 


is  equivalent  to  the  computation 


iwkm< 


+  a2,  2 

Y  X 


/  2X  2.  2 
wCY  ‘  m  X 
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C3Y 

2  ”7,2  ? 
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msXSY 


for  i  =  1,2,  *  is  the  PtZ_<  z]  when  2  is  a  standard  normal 

random  variable,  and  ?  is  the  density  function  of  a  standard 
normal  random  variable. 

If  the  substitution  Y  =  aX-c  were  made  in  the  density 

f„  v(x,y),  an  analogous  situation  would  arise  with  the  limits 
"  *  * 

of  the  integration  being  X.  .  and  X.  _ . 

7.  If  the  value  of  the  density  under  the  bivariate 
normal  curve  over  the  region  of  the  ellipse  defined  by  the 
subinterval  £5^  is  weighted  by  the  probability  of  observing 
a  bearing  error  (5^.-5)  the  result  is  the  probability  of 
observing  the  bearing  of  the  threat  as  3  when  the  true  threat 
location  is  in  the  segment  of  the  ellipse  defined  by  ££,^. 

This  probability  is  approximated  as  follows: 


.1 I!"  VPqlMWM  . . . .  WWff, 


?k  '  fX,Y(^i“'y)  dy  I  fg{3;6k)* 

Ykl 


Since  the  bearing  error  is  assumed  normally  distributed  with 

2 

mean  zero  and  variance  ,  the  value  of  fg(g;S^)  can  be 


determined  by  the  expression 


1  8k  "  6 
~~  t  ( JS- - ) , 

a  8  °8 


where  <p  is  as  above  the  density  function  of  a  standard 
normal  random  variable. 


f 1 (3)  =  /  fft(S;8  )  fa  (3  )  dB 

all  3,  B  PmorTTT 


'TRUE 


is  approximated  by 


2n  *  j^2 

l  IV  /  fa(8.-3t) 


L  ■  V  1 

k=l  K  y 


'sl*'uk' 


The  value  of  the  sum  is  determined  by  repeating  steps 
A. 6  and  A. 7  above  at  the  midpoints  of  each  of  the  2n  sub¬ 
intervals  defined  in  step  A. 4  and  summing  the  results  of 
each  of  these  calculations.  Obviously,  the  finer  the  sub- 
divisicu  of  (3^,3^),  the  more  accurate  will  be  the  estimate 
of  the  likelihood,  but  also  the  longer  the  calculation  will 
take. 


B.  SENSOR  POSITION  UNCERTAINTY 

The  result  of  the  above  calculation  will  be  the  likeli¬ 
hood  f'(3)  that  the  threat  lies  along  the  bearing  measured 
by  the  sensor  given  that  the  threat  is  located  within  the 
threat  ellipse  and  the  true  bearing  is  within  the  interval 
(8T,6r.)  and  the  sensor  is  at  the  position  used  to  perform 

i-i  u 

the  calculation.  We  next  will  introduce  additional  calcu¬ 
lations  that  are  required  to  account  for  the  fact  that  the 
sensor  position  is  not  known  with  certainty. 

1.  If  the  position  of  the  sensor  is  estimated  as  being 

within  R  nautical  miles  of  (uo,vQ),  its  assumed  coordinates 

in  the  xy-system,  with  p3  *  100%  confidence,  estimate  the 

parameters  of  the  sensor  location  density  fy  V(u,v)  with 

2 

mean  zero  and  variance  a  :  mean  =  (uo,vQ) ,  variance, 

2  2 

a  -  R  /  (-2  ln(l-p3))  in  threat  centered  coordinates. 

2.  The  bearing  measured  by  the  sensor  is  S  regardless 

of  the  sensor  location.  Assume  that  the  area  of  intersection 
of  the  angular  wedge  (e^Sy)  and  the  threat  ellipse  does  not 
change  significantly  as  the  sensor  position  is  moved  along 
the  line  of  bearing  3.  Then  the  effect  of  the  bivariate 
normal  distribution  of  sensor  location  can  be  approximated 
by  considering  only  the  univariate  normal  density  along  a 
line  through  (uo,vQ)  perpendicular  to  3.  Repeat  the  calcu¬ 
lations  in  steps  A. 3  through  A. 8  above  with  the  sensor  located 
at  each  of  the  three  points  {uo,vQ) ,  (uo+.970COs(9-S-9O) , 

v  +. 97osin (e-0-90) ) ,  and  (u  - .S7ocos (e-S-90) ,v  -97asin(6-S-90) ) . 
o  00 

If  the  line  through  (uQ,v0)  perpendicular  to  0  is  subdivided 
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symmetrically  about  (uo,vo)  such  that  1/3  of  the  univariate 
normal  density  lies  above  each  subinterval,  the  three  points 
chosen  above  represent  the  "center  of  gravity"  of  each  third 
of  the  density  (Figure  6) . 

3.  If  R  is  chosen  to  include  a  significant  proportion 
of  the  sensor  density,  i.e.,  on  the  order  of  2c  or  greater, 
the  probability  of  the  sensor  being  located  in  each  of  the 
three  regions  is  approximately  1/3.  Thus,  if  p^  is  on  the 
order  of  .86,  multiply  each  result  in  step  B.2  by  1/3. 

4.  Summing  the  results  of  steps  B.2  and  B.3  yields  an 
estimate  of  the  likelihood  that  the  threat  lies  at  bearing 
3  given  the  threat  is  in  the  p^  *  100%  confidence  ellipse, 

8t  is  in  (3^,8^)  and  the  sensor  is  in  a  p^ *  100%  confidence 
region.  That  is,  the  likelihood  is  estimated  by 

3  i 

f'{3)  =  l  i  f,'(3)  , 

j=l  J 

i 

where  j  is  the  index  of  sensor  position  in  figures  5  and  6. 

Instructions  for  application  of  the  TI-59  program  to 
calculate  the  above  are  included  in  Appendix  A. 

C.  THREAT  DISTRIBUTION  NOT  BIVARIATE  NORMAL  AFTER  TIME 

LATE  ELAPSED 

The  basic  approach  to  evaluating  f'(8)  when  the  time  late 
distribution  of  the  threat  is  no  longer  bivariate  normal 
is  the  same  as  that  just  discussed.  The  principal  difference 
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FIGURE  6.  ESTIMATE  OF  SENSOR  POSITION  DENSITY 


arises  from  the  fact  that  the  jne  late  density  is  signi¬ 
ficantly  more  complex  in  this  case. 

The  threat  distribution  becomes  distorted  from  the  normal 
after  some  time  late  when  the  course  and  speed  of  the  threat 
are  constant  but  not  known  with  certainty.  Application  of  the 
method  described  herein  requires  that  the  user  assume  a  dis¬ 
crete  distribution  of  the  speed  of  the  threat  with  upper  and 
lower  bounds  SM  and  respectively.  In  addition,  the  threat 
course  is  assumed  to  be  uniformly  between  0°  and  360°. 

The  density  of  threat  location  after  some  time  late  t^ 
when  the  threat  speed  is  s^  then  becomes  [Ref.  1] : 


.  360  .  x-s. t^cos(e-ij;) ) 

f  (x,y;t  ,s.)  =  v-~~  --  /  exp[-5-( - 5 - 

i  •  cx°y  0  •* 


X 


ly-siVin(e-*))  ^  d„ 

+  _  > J  36Q 


where  9  is  the  bearing  of  the  major  axis  from  North.  Note 
that  the  new  threat  density  is  still  centered  at  the  same 
position  as  the  time  zero  ellipse  but  its  shape  changes  as 
in  Figure  4  of  Chapter  II.  If  the  threat  speed  is  s.,  and 
the  course  is  uniformly  distributed  over  (0o,3606),  the 
outer  limit  of  the  new  planar  region  containing  the  threat 
after  time  has  elapsed,  given  that  it  was  originally 
located  in  the  p^  *  100%  confidence  ellipse  with  semi-axes  A 
and  B,  can  be  represented  by  an  ellipse  with  semi-major  axis 
A+s . t_  and  semi-minor  axis  B+s.t_  (Figure  7).  Thus  the 

1  L  X  Lj 
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PLANAR  REGION 


region  over  which  the  density  will  be  evaluated  is  still 
the  intersection  of  an  ellipse  with  an  angular  wedge, 
recalling  however  that  the  density  is  no  longer  normal. 
Further,  the  new  elliptical  region  does  not  represent  a 
p^  x  100%  confidence  region  of  the  time  late  density. 

Calculation  of  the  likelihood  that  the  threat  lies  along 
bearing  6  given  that  the  threat  lies  in  the  p^  *  100%  ellipse 
at  time  zero,  that  the  true  bearing  lies  in  (3^,5^),  that 
the  sensor  lies  in  the  p^  *  100%  circle  and  that  the  speed  is 
s^  proceed*  as  follows: 

1.  Estimate  the  parameters  of  the  original  normal  dis- 

2  2 

tribution:  mean  =  (0,0),  =  A  / (—2  ln(l-p^)), 

oy2  =  B2/(-2  ln(l  -Pl))  . 

2.  Determine  the  upper  and  lower  bounds  on  the  true 
bearing  wedge,  3  +  kgOg  and  8-kaoa. 

3.  Determine  the  time  late  planar  region  as  the  ellipse 
with  semi -major  axis  equal  to  A+s.tT  and  semi -minor  axis 
equal  to  B  +  s^t^. 

4.  Determine  the  position  of  the  sensor  relative  to  the 
ellipse  center. 

5.  Subdivide  the  bearing  fan  (3-knai,,3  +  kac0)  into 

o  e  p  3 

2n  subintervals. 

6.  At  the  midpoint  of  each  subinterval,  determine  the 
equation  of  the  line  through  the  ship  position  at  that 
bearing  S^. 

7.  Let  the  equation  of  the  line  of  bearing  3^  be 
Y  -c 

X  =  — — .  Then,  the  plane  perpendicular  to  the  xy-plane  which 


Wv  and  the  area  under  the  curve  gT  between  Y, ,  and  Y, - 


The 


area  under  curve  gL  can  be  approximated  as  follows.  Sub¬ 
divide  the  interval  (Yj^Y^)  ^nto  n-s  segments  of  length  h. 


Evaluate  gT  at  the  midpoint  of  each  segment,  y_.  Note  that 

Ii  U 


the  first  term  in  parenthesis  in  (4)  is  a  constant  equal  to 


m  times  the  value  of  the  densitv  of  a  normal  random  variable 


with  mean  zero  and  variance  +  evaluated  at  c.  The 


second  term  in  parenthesis  is  the  density  function  of  a  normal 
random  variable  with  mean 


and  variance 


CCY2 
+  tv?  a 
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2  2  2 

m  °X  °Y 
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evaluated  at  y.  The  variable  y  also  appears  in  the  integral 
term  in  equation  (4).  Numerically  evaluate  this  term  of 
(4)  with  y  =  y  .  Let  gT  (y  )  be  the  result  of  this  computation 
Then,  evaluating  the  area  under  the  curve  gL  between  Ykl 
and  Yk2  is  equivalent  to  the  calculation: 
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v3y  +m  “X 


£ 


-) 


/ 


2  2  2 
ay  +  ra  ox 


2  2  2 
Y  +m  °X 


mo  ^Oy 


gI(yJ) 


where  ${*)  is  the  density  function  of  a  standard  normal  random 
variable,  and 


J 


l  ^  £  £ 

5y  +  m  °x 


mc^Oy 


An  analagous  situation  arises  if  the  substitution  Y  =  mX  +  c 
is  made  for  y  in  the  time  late  density  fx  ^{X/yyt^} . 

8.  The  value  of  the  density  over  the  subinterval  con¬ 
taining  3^  is  then  weighted  by  the  instantaneous  probability 
that  the  bearing  error  is  3^  -  3: 


Wk  v  J  fX,Y(XL'yL;tX)dyLi  fS(8?Sk) 
Xkl 


9.  Steps  C.  7  and  C.8  are  then  repeated  for  each  sub¬ 
interval  of  (s  ,3„)  and  the  results  of  each  calculation  are 
Xj  u 

summed. 

10.  The  uncertainty  in  sensor  location  is  accounted  for 
by  repeating  steps  C.4  through  C.9  with  the  assumption  the 
sensor  is  located  at  each  of  the  three  points  in  Figure  5, 
multiplying  by  the  probability  the  sensor  lies  in  that 
interval  and  summing  each  result. 

11.  The  result  of  calculations  in  steps  C.l  through 
C.10  is  the  likelihood  that  the  threat  lies  on  bearing  3 


given  the  threat  was  originally  located  in  the  p^  x  100% 
confidence  ellipse,  3^  z  (3^,3^},  the  sensor  is  located  in 
the  p^  x  100%  circle  and  the  threat  speed  is  s^.  The  con¬ 
dition  that  the  soeed  is  s .  is  removed  by  repeating  the 


calculations  C.l  through  C.10  for  each  of  the  speeds  s^, 
i  =  1, . . . ,M  multiplying  the  result  by  the  probability  that 
the  speed  equals  s^  and  summing  all  M  results.  The  final 
result  is  the  likelihood  that  the  threat  lies  on  3  given 
it  was  originally  located  in  the  p^  *  100%  ellipse, 

3j,  e  and  the  sensor  is  in  the  p^  *  100%  circle  esti¬ 

mated  as 


[  P[S=s.]  I  \  f  •  (3)  . 
i=l  *  j=l  J  3 
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IV.  CONCLUSIONS 

Possible  applications  of  and  extensions  to  the  algorithm 
developed  in  Chapter  III  are  discussed  in  this  chapter. 

\s  indicated  in  Chapter  I,  the  objective  of  this  paper 
has  been  to  develop  a  means  of  assessing  the  likelihood  a 
threat  whose  position  has  previously  been  estimated  lies  at 
a  given  measured  bearing  from  a  local  sensor.  The  procedures 
were  developed  with  a  view  towards  permitting  the  user  to 
make  comparisons  among  lines  of  bearing  measured  by  different 
sensors  or  among  conflicting  bearing  information  generated 
by  one  sensor.  The  approach  chosen  has  been  to  estimate  the 
likelihood  that  a  threat  lies  on  bearing  a  given  that  the 
threat  is  located  in  an  ellipse  of  specified  confidence 
P1  x  100%,  that  S  is  measured  with  p2  *  100%  accuracy,  and  that 
the  sensor  position  is  measured  with  p^  *  100%  accuracy.  The 
algorithm  to  calculate  the  likelihood  in  the  cases  where  the 
probability  distribution  of  the  target  can  be  assumed  to  be 
bivariate  normal  at  the  time  of  the  bearing  measurement  has 
been  implemented  on  the  TI-59  calculator.  The  cases  in  which 
this  program  applies  sure  the  following:  (1)  when  the  time 
elapsed  since  generation  of  the  threat  error  ellipse  is  small 
enough  to  justify  using  the  original  estimate  of  the  ellipse; 
(2)  when  threat  course  and  speed  are  known,  in  which  case 
the  ellipse  center  is  translated  from  the  original  position 
according  to  the  course,  speed  and  elapsed  time  information; 
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(3)  when  the  threat  can  be  assumed  to  be  moving  about  in  a 
random  manner  over  a  significant  region  in  such  a  way  that 
the  ellipse  center  remains  unchanged,  but  the  x  and  y  vari¬ 
ances  have  increased.  When  none  of  these  cases  hold,  but 
the  course  is  assumed  uniformly  distributed  over  (0°,350°) 
and  speed  has  a  discrete  distribution  over  a  finite  interval, 
the  threat  density  at  the  time  late  t^  is  not  a  bivariate 
normal  density.  The  algorithm  applicable  in  this  case  has 
not  been  implemented,  but  is  described  in  some  detail  in 
Chapter  III  and  Appendix  B. 

Once  the  appropriate  computation  has  been  completed  for 
each  of  the  bearings  considered,  the  results  can  be  used  to 
weight  the  value  of  several  bearings  in  refining  the  threat 
location  estimate  provided  by  the  external  sensor.  Note  that, 
although  unlikely,  f ' { 3)  may  correctly  be  greater  than  one. 
Comparisons  using  these  likelihoods  should  be  made  only  when 
the  upper  and  lower  bounds  on  the  true  bearing  fan  for  each 
bearing  are  chosen  at  the  same  probability  level  p2  and  the 
uncertainty  areas  for  all  sensors  include  the  same  probability 
level  p^.  Further,  these  likelihood  levels  should  be  selected 
so  as  not  to  exclude  a  significant  portion  of  the  appropriate 
density.  If  p2  or  p^  are  not  the  same  in  all  cases  to  ue 
compared,  f 1 (5)  must  be  divided  by  the  applicable  value  of 
P2  or  p.j  for  each  measured  bearing  3  to  be  considered. 

Having  established  the  relative  value  of  available 
bearing  information,  the  user  can  allocate  weapons  or  further 
search  effort  accordingly.  However,  the  probabilities  calculated 
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are  strictly  ordinal  data  and  do  not  define  a  redistribution 
of  target  location  probability  based  on  additional  informa¬ 
tion.  Further,  the  threat  ellipse  does  not  contain  the 
target  with  certainty.  The  power  to  predict  the  probability 
of  success  of  a  search  or  weapons  allocation  plan  based  on 
the  priorities  established  by  these  procedures  is  limited 
by  these  constraints.  In  this  area  in  particular  further 
research  would  be  useful. 

In  situations  such  as  that  for  which  the  procedures  in 
this  paper  have  been  developed,  where  a  track  has  not  been 
developed  on  the  target,  introduction  of  unrelated  bearing 
data  to  a  target  motion  model  could  impact  significantly  on 
Mie  reliability  of  future  position  predictions.  If  there  is 
high  confidence  in  the  reliability  of  the  estimate  of  the 
threat  ellipse  provided  by  the  ocean  surveillance  sensor, 
the  prioritization  established  herein  could  be  used  to  pro¬ 
cess  bearing  data  prior  to  input  to  a  target  motion  model. 

Using  a  pre-established  threshold,  only  those  bearings  which 
coincide  with  the  threat  with  an  acceptable  level  of  like¬ 
lihood  could  be  used  to  refine  or  update  a  track  on  the 
threat . 

Desirable  enhancements  to  the  algorithms  include  providing 
for  the  instances  in  which  the  interval  of  uncertainly  of 
the  target  course  is  known  to  be  less  than  (0s, 360°).  Further, 
if  the  circular  region  of  radius  R  contains  the  sensor  with 
significantly  less  than  86%  confidence  or  the  assumption  of  a 
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bivariate  normal  distribution  of  sensor  location  is  unsatis¬ 
factory,  it  is  left  to  the  user  to  modify  the  calculations 
accordingly. 

The  utility  of  these  algorithms  would  also  be  improved 
by  implementation  on  a  larger  and  faster  system  than  the 
TI-59  calculator. 
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APPENDIX  A.  TI-59  PROGRAM  VERBAL  FLOW  AND  USER'S  INSTRUCTIONS 


Part  I 
Step  Number 

000  -  029 

020  -  029 

030  -  052 

053  -  058 

059-063 

064  -  073 


Verbal  Flow 

Enter  the  confidence  level  for 
the  threat  ellipse.  Calculate  the 
value  of  k  for  the  given  p^: 

k  =  /-2  ln(l-p^)  . 

Enter  the  length  of  the  semi-major 
axis,  A.  Calculate  ax  =  A/k. 

Enter  the  length  of  the  semi-major 

axis,  B.  Calculate  o  =  B/k. 

*  2  2 

Calculate  B/A,  axoy,  ax  and  oy  . 

Enter  orientation  of  semi-major 
axis,  d. 

Enter  bearing  from  sensor  to  center 
of  threat  ellipse,  a.  Calculate 
9— a . 

Enter  distance  r  from  sensor 
position  to  center  of  threat 
ellipse.  Determine  rectangular 
coordinates  of  sensor  position 
( U, V)  from  polar  coordinates 
( -r , e-a)  . 


074  -  089 


Store  the  constants  360,  /2-rr  . 
Initialize  register  35  to  0. 


Enter  number  of  standard  devia¬ 
tions  desired  for  bearing  fan,  k. . 

6 

Enter  standard  deviation  of  bearing 

error,  c,.  Calculate  k  o  . 

3  S3 

Enter  the  angular  stepsize  desired, 
AS ,  for  incrementally  stepping 
through  (Sj^Sy)  .  Calculate  jj-AB  . 
Calculate  the  largest  number  n  of 
increments  of  size  AS  contained  in 
k3a.  “  degrees.  Initialize 

counter  00  to  n+1.  Save  n+1  in 
register  20.  Determine 

5 AS  =*  k„oa  -  iAS  -  nA3  , 

-  b  b  *■ 

the  residual  increment.  Calculate 

j<$A3  . 

Initialize  counter  01  to  2. 

Calculations  will  be  made  at  the 

midpoint  of  each  interval  from 

S  +  4as  to  S  +  k3o,,  then  at  the 
i  3  b 

midpoint  of  each  interval  from 
3  -  jAS  to  3  -  k^o^  »  and  finally 
at  3.  Counter  01  indicates  whether 
calculations  are  complete  on  both 
sides  of  6. 

Enter  bearing  measured  by  sensor, 

3. 

Enter  standard  deviation  of  sensor 
position  a. 


2 


155 


Enter  index  of  the  sensor  position 
to  be  used  for  this  run. 


156  -  157 


158  -  238 


239  -  259 

260  -  340 


Coordinates  of  the  sensor  position 
are  selected  in  accordance  with 
run  number  entered  above  in 
Subroutine  sin. 

Determine  whether  a  bearing  S ' 
parallel  to  either  axis  is  included 
in  the  fan  (SL,Sy).  If  (SL,Sy) 
includes  a  bearing  parallel  to  the 
y-axis,  use  program  1  for  Part  II. 
The  appropriate  program  number  is 
displayed  in  calculator  display 
register.  If  (Sy,3y)  does  not 
contain  a  bearing  parallel  to 
either  axis,  use  program  1. 

Subroutine  P/R. 

Subroutine  sin. 
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Part  II,  Programs  1  and  2 

Step  Number  Verbal  Flow 


000  -  003 


Initiate  3 ' 


004  -  007 


If  the  last  angular  increment  on 
this  side  of  S  has  been  considered, 
go  to  step  513-  Otherwise 
continue. 


008  -  012 


Decrement  counter  for  angular 
increment.  If  counter  =  0,  go  to 
021.  Otherwise  continue. 


013  -  020 


021  -  028 


Remove  flag  to  indicate  this  is 
not  the  last  angular  increment. 
Recall  43,  the  input  angular  step- 
size.  Go  to  029. 

Set  flag  to  indicate  this  is  the 
last  angular  increment  on  this 
side.  Add  one-half  the  input 
angular  stepsize  43  and  one-half 
the  residual  stepsize  543: 

^■48  +  ^-5  43. 


029  -  030 


Increment  B '  by  the  appropriate 
stepsize. 


031  -  O'* 5 


Convert  S’  to  an  angle  between  O' 
and  360°. 


046  -  051 


Calculate  9-8’.  Print  9-3'. 
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If  ie-S'i  =  90°  or  =  270°,  go  to  236 


If  ! e-S '  |  =  0°  or  =  180°,  go  to  075 
Otherwise  go  to  133. 

If  the  absolute  value  of  the  y 
coordinate  of  sensor  position  is 
greater  than  the  length  of  the 
semi-minor  axis,  £,  go  to  004.  In 
this  case  3 '  does  not  intersect  the 
error  ellipse.  Otherwise  continue. 

Set  flag  2  to  indicate  that  the 
bearing  3 '  is  parallel  to  the 
x-axis . 

Calculate  the  coordinates  of  the 
points  of  intersection  of  3'  with 
the  threat  ellipse.  The  y-coordi- 
nates  are  equal  to  V,  the 
y-coordinate  of  sensor  position, 
x-coordinates  are  determined  in 
Subroutine  y  .  The  points  of 
intersection  are  symmetric  about 
the  y-axis.  Thus,  =  -X^2- 
Store  the  smaller  x  value  in 
register  27,  the  larger  in  register 
28. 

Save  locations  of  X^,  Xk2,  Y^, 
ax,  and  c Y-  Go  to  291. 

Remove  flags  2  and  3  to  indicate 
that  3 '  is  not  parallel  to  either 
the  x-  or  y-  axes . 


Calculate  the  slope  of  3 ' : 
m  =  tan  (9-3  ' )  . 

Calculate  the  y-intercept  of  3 ' : 
c  =  V  -  mU. 

If  >  A^m^  +  B^,  go  to  004.  In 
this  case,  3*  does  not  intersect 
the  error  ellipse.  Otherwise 
continue. 

Calculate  Xkl  and  Xk2,  the  x- 
coordinates  of  the  points  of  inter 
section  of  3*  with  the  ellipse: 


Calculate  the  y-coordinates  of  the 
points  of  intersection  of  3'  with 
the  ellipse,  Ykl,  Yk2: 


Yki  =  "*ki  +  c 


Yk2  mXk2  +  c  ' 


Go  to  299. 


If  luj,  the  absolute  value  of  the 
x-coordinate  of  sensor  position, 
is  greater  than  the  length  of  the 
semi-major  axis.  A,  go  to  004.  s' 


does  not  intersect  the  error 
ellipse  in  this  case.  Otherwise 
continue. 


248  *  249 


250  -  272 


273  -  298 


299  -  328 


Set  flag  3  to  indicate  that  3  *  is 
parallel  to  the  y-axis. 

Calculate  the  coordinates  of  the 
points  of  intersection  of  3 '  with 
the  error  ellipse.  The  x- 
coordinates  are  equal  to  U,  the 
x-coordinate  of  sensor  position. 
The  y- coordinates  are  determined 
in  Subroutine  yX.  The  points  of 
intersection  are  symmetric  about 
the  x-axis.  Thus,  Y,  .  =  -Y,  _. 
Store  the  smaller  y  value  in 
register  29 ,  the  larger  in  30 . 


Save  the  locations  of  Y^2,  Y^, 


Tel1 


X* 


Determine  the  width  of  the  strip 
around  S',  W: 


If  flag  0  set. 


,1. 


W  =  (d^+d2)  tanf-joAB), 


Otherwise 


W  =  (d1+d2)  tan  ^AB) . 

d-^  and  d2  are  the  distances  from 
(U, V) ,  the  sensor  position,  to  the 
intersection  points  (X^i  ,Y^)  and 
(Xk2,Yk2} "  respectively.  Distances 
are  calculated  in  Subroutine  log. 
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Multiply  W  by  S3  if  flag  0  not 
set.  Otherwise,  multiply  by  5 AS. 
Save  result  in  register  37. 

329  -  334  If  S'  is  parallel  to  either  the 

x-  or  the  y-axis,  go  to  468. 
Otherwise  continue. 


When  using  program  1: 


335  -  380 


381  -  391 


392  -  401 


Express  the  v-coordinates  of  the 
intersection  points  as  standard 
normal  random  variables,  Z^  and  Z2 


Yki  - 


Cay 

„  2  2  _  ; 
*  Y  +  m  '"'X 


iiU  , 


XaY 


r  ->  22 

T /cv‘  +  m  ~ 


i  =  1,2. 


Sort  the  values  Z^  in  descending 
order.  Let  Z^'  be  the  larger 
value.  Z^'  is  the  smaller. 

Calculate  ( Z2 ' )  ~  the 

probability  that  a  standard  normal 
random  variate  lies  between  Z^' 
and  Z2'. 
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402  -  424 


Multiply  the  results  of  steps 
392-401  by: 


The  result  o 


When  using 

335  -  372 


373  -  383 


/  2  2  2  ,/  2“  2“  2 
V - y  +  m  cx  i-cv  ^  m  ~X 


-}x:  ‘  S3  ’  x  W, 


where  :  is  the  standard  normal 
density  t motion. 


f  the  calculations  in  steps  335-424  is 


y:<2 


iw 


y 


kl 


rogram  2: 


Express  the  x-coordinates  of  the 
intersection  points  as  standard 
normal  random  variables  Z-^  and  Z2- 


Z. 

l 


me; 


Xk2  + 


X 


/  T~  2  2 

1/c  -r  m  c? 


x-v 


./  2  2  2 
1/  +  m  a 


I 


i  =  1,2. 


Sort  the  values  of  Z^  in  descending 
order.  Let  Z2 '  =  larger  value. 

Z^'  =  smaller. 


384  -  393 


Calculate  MZ2’)  -  S(Z^')i  the 
probability  that  a  standard  normal 
random  variate  lies  between  Z^' 
and  Z2 ' . 

3?4  -  413  Multiply  the  results  of  steps 

384-393  by: 


{- 


2  .  2  2 
v  "X 


l)cv2  +  m2c2 


-)  *  W, 


where  ?  is  the  standard  normal 
density  function. 


The  result  of  the  calculations  in  steps  335-413  is 

Y, 


k2 


jW 


fy.  Y(x,mx+c)  dxj 


Regardless  of  which  program  is  in  use: 


425  -  465 


Multiply  results  of  previous 
calculation  by 


1  .  .3' 

-—i'k— 


where  S' 


2  ,  B 

3  * -360 0  ,  S'  >  180° 


and  ;  is  the  standard  normal  density. 


The  result  of  this  calculation  is 


f,(5;3'}  jW  j  fx  y(x,y)  dyj 
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466  -  467 


Go  to  508 


468  -  507  Calculate  the  area  under  the  curve 

formed  by  the  intersection  of  = ’ 

with  fY  v(x,y).  When  S'  parallels 

a,  i. 

the  x-axis  this  calculation 
becomes : 

—  s(— )  [i  (Z  )  -  i(Z  )  ]  , 

~Y  “Y  z  x 

where  Z2  =  \2^X  and  Z1  =  ^Scl^'X  ' 
V  is  the  y-coordinate  of  the  sensor 

position. 

When  s  *  is  parallel  to  the  y-axis , 
the  calculation  is: 

[i(Z_)  -  c (Z, )  ]  , 

■'X  2  1 

where  Z-  =  Yk2/;y  and  =  Ykl/cY  , 
U  is  the  x-coordinate  of  the  sensor 
position. 

Go  to  422  to  complete  calculation 
of  f. (£;=’) !W  / fx  y{x.y) i  . 

(Go  to  410  in  Program  2) . 

508-512  Accumulate  the  probability  at  each 

angular  interval.  Display  result. 

513  -  515  If  flag  1  is  set  go  to  558. 

Otherwise  continue. 

516  -  520  If  flag  0  is  not  set,  that  is  if 

the  angular  increment  just  con¬ 
sidered  was  not  the  last  on  this 
side  of  5,  go  to  004  and  continue 
»  calculation  on  same  side  of  3. 

Otherwise  continue. 


51 


521 


524 


Decrement  the  counter  01.  If 
counter  is  now  0  go  to  543.  In 
this  case,  the  probabilities  have 
been  calcic- .ed  on  both  sides  of  £. 
The  calculation  at  s  remains  to  be 
done.  Otherwise  continue. 

525  -  542  Remove  flag  0.  Multiply  is,  ^ls , 

and  jfiS  by  -1.  Reinitialize 
counter  00  to  n+1.  Go  to  000  to 
begin  calculation  on  second  side 
of  5 . 


543  -  567  Set  flag  01  indicating  that  calcu¬ 

lations  on  both  sides  of  £  have 
been  completed.  The  next  iteration 
will  do  the  calculation  at  5  *  =  5 . 
Remove  flag  00.  Initialize  £  *  to  0 
Recall  £ .  Go  to  029. 

558  -  566  Remove  flag  01.  Display  the 

accumulated  likelihood.  STOP 


The  result  of  this  calculation  is  approximately  £_,'{£ } 

J 

sensor  at  (U. ,V. }  where  j  =  run  number  entered  in  Part 


wrtn 


567  -  587  Subroutine  yx  calculates  the  points 

of  intersection  when  s  *  is  parallel 
to  either  the  x  or  v  axis: 


Parallel  to  y-axis: 


58S  -  610  Subroutine  log  calculates  the 

distance  from  sensor  position  to 
point  of  intersection  of  S ’  with 
the  ellipse: 
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part  III 
Step  Number 


Verbal  Flow 


000  -  012 


013  -  021 


022  -  029 


030  -  032 


Enter  result  of  run  1. 
Multiply  by  1/3. 

Enter  result  of  run  2. 
Multiply  by  1/3. 

Enter  result  of  run  3. 
Multiply  by  1/3. 

Display  likelihood. 


033-041  If  p2  is  the  same  for  all  bearings 

to  be  compared,  enter  1.  Go  to 
042.  Otherwise,  enter  p2  for  this 
bearing.  Divide  likelihood  by  p2- 
Display  result. 

042  -  050  If  p^  is  the  same  for  all  bearings 

to  be  compared,  enter  1  and  STOP. 
Otherwise,  enter  p3  for  this  bearing. 
Divide  likelihood  by  p^.  Display 
result.  STOP. 


USER'S  INSTRUCTIONS 

The  program  to  determine  the  likelihood  that  the  threat 
lies  along  bearing  S  given  the  threat  is  in  the  confidence 

3  x  100% 

confidence  circle  is  in  three  parts.  All  parts  require  the 
use  of  a  printer  and  the  use  of  the  Applied  Statistics 
Library  Module.  Prior  to  running,  the  calculator  must  be 
repartitioned : 

1.  Enter  4 

2.  Press  2nd  OP  17 

Part  I 

1.  Read  sides  1  and  2  of  Part  I 

2.  Read  side  4  of  Part  II  either  program  1  or  program 
2.  Since  program  2  is  used  more  often,  unless,  it  is  known 
that  the  bearing  fan  (3^3^,)  contains  a  bearing  parallel  to 
the  major  axis,  recommend  using  program  2  of  Part  II. 

3.  Press  R. 

4.  Enter  p^,  the  confidence  level  of  the  threat  ellipse. 
Press  A 

5.  Enter  A,  length  of  the  semi-major  axis.  If  data 
provided  is  length  of  the  entire  major  axis,  divide  by  2 
before  entering. 

Press  B. 

6.  Enter  B,  length  of  the  semi-minor  axis.  If  data 
provided  is  length  of  entire  minor  axis,  divide  by  2  before 
entering.  Press  C. 


ellipse,  ^  s  (B^Sy)  and  the  sensor  is  within  the  p 


7.  Enter  9,  the  bearing  of  major  axis  from  North. 

0°  <  0  <  180°.  9  is  entered  in  degrees.  Press  D. 

8.  Enter  a,  the  bearing  of  the  threat  ellipse  center 
from  the  estimated  sensor  position.  0°  <_  a  <_  360°,  in 
degrees.  Press  E. 

9.  Enter  r,  the  range  from  the  estimated  sensor  posi¬ 
tion  to  the  center  of  the  threat  ellipse.  Press  2nd  A'. 

10.  Enter  k  ,  the  number  of  standard  deviations  desired 

3 

in  one  direction  from  3.  The  program  will  construct  a  fan 
of  equal  size  on  the  other  side  of  3.  Press  2nd  B'. 

11.  Enter  a.,  the  standard  deviation  of  bearing  error. 

B 

Press  2nd  C' . 

12.  Enter  AS,  the  desired  angular  stepsize.  Press  2nd 

D' . 

13.  Enter  s,  the  measured  bearing.  Press  2nd  E*. 

14.  Enter  a,  the  standard  deviation  of  sensor  position. 
Press  R/S. 

15.  Enter  the  number  of  this  run,  1,  2  or  3.  When  run 

number  =  1,  sensor  is  located  at  its  estimated  position 

(uQ,v0) .  When  run  number  =  2,  the  location  will  be 

(u  + .  97o*cos  (0-3-90)  ,v +.  97asin  ( e-s-90)  )  .  When  run  number  = 
o  o 

the  location  will  be  . 97acos  (e-s-90)  ,v0-.  970sin  (6-6-90) ) 
Press  R/S.  If  program  number  displayed  matches  that  of  the 
Part  II  side  4  read  in,  continue  to  16.  Otherwise,  press 
2nd  CMS,  RST.  Read  in  side  4  of  the  Part  II  program  which 
corresponds  to  the  number  displayed.  Repeat  3  through  15. 


16.  Press  4  2nd  WRITE.  Rerecord  side  4  of  the  Part  II 
program  read  in.  This  enables  data  entered  in  Part  I  to 
be  transferred  to  Part  II. 

Part  II 

1.  Read  sides  1,  2,  3  and  4  of  Part  II  Program  1  or  2 
as  selected  by  the  Part  I  program. 

2.  Press  RST 

3.  Press  R/S.  Values  of  3*.  and  Wjf  (S;S.)f  (S.) 

v.  S  t  STRUE  t 

will  print  alternately.  Final  result  also  prints  out  at  end 
of  calculation. 

4.  Record  final  result:  Likelihood  threat  lies  along 
bearing  3  given  threat  is  in  p^  *  100%  ellipse,  true  bearing 

is  in  p2  x  100%  fan  and  sensor  is  at  location  used  in  this  run. 
Parts  I  and  II  must  be  completed  3  times  (Run  numbers  1,  2 
and  3)  before  proceeding  to  Part  III  if  uncertainty  in  sensor 
position  is  being  considered. 

Part  III 

1.  Read  side  1  of  Part  I.  I. 

2.  Enter  result  of  run  1.  Press  A. 

3.  Enter  result  of  run  2.  Press  B. 

4.  Enter  result  of  run  3.  Press  C. 

5.  Enter  p2,  confidence  level  of  (S^Sy),  if  necessary. 

Otherwise  enter  1.  Press  D. 

6.  Enter  p^  confidence  level  of  sensor  position,  if 
necessary.  Otherwise,  enter  1.  Press  E. 


APPENDIX  B.  TIME  LATE  VERBAL  FLOW 


VERBAL  FLOW 

1.  Enter  p^.  Determine  the  value  of  k:  k  =  /-2  In  (1  -  p-^) 

2  2 

2.  Enter  A.  Determine  ox  =  (A/k) 

2  2 

3.  Enter  B.  Determine  =  (B/k) 

4.  Enter  e,  orientation  of  major  axis 

5.  Enter  a,  bearing  from  sensor  to  ellipse  center 

6.  Enter  range  r  from  sensor  to  ellipse  center 

7.  Calculate  (u  ,v  )  the  coordinates  of  mean  sensor 

o  o 

position.  uQ=*  -r  cos(e-o),  vQ=  -r  sin(e-a) 

8 .  Enter  a . ,  the  standard  deviation  of  bearing  error 

D 

9.  Enter  k  ,  the  number  of  standard  deviations  to  be 

3 

included  in  the  bearing  fan  on  each  side  of  3. 

10.  Enter  desired  angular  stepsize,  A3 

11.  Calculate  the  number  of  iterations  of  size  A3  required 

on  each  side  of  3:  I  =  [ (k.o -  iA3)/A3]  where  [•]  means 

00  *• 

the  greatest  integer  less  than  or  equal  to  the  value  within. 

12.  In  oeneral  (k.o.  -  4as)/a3  is  not  an  integer.  Determine 
the  size  of  the  fractional  increment: 

5 AS  =  ( (kgag  -  |as)/A3  -  I) AS. 

13.  Enter  3/  the  measured  bearing.  Let  PROB2=0,  PROB1=0 

14.  Enter  the  time  late  t^ 

15.  For  each  of  the  discrete  threat  speeds  to  be  con¬ 
sidered,  repeat  steps  16  to  43.  Then  go  to  44. 

16.  Enter  the  target  speed  s.  Let  PROB  =  0 


17.  Expand  the  outer  limit  of  the  threat  ellipse: 

Let  A*  =  A  +  st^ 

B’  =  B  +  stL 

18.  Let  S'  =  3 

19.  For  each  of  the  I  increments  of  size  A3  repeat 
steps  20  to  38. 

20.  Let  S'  =  S'  +  AS 

21.  Calculate  (3-3’).  If  ■  e-3 '  j  *  0°  or  180°,  go  to  27 
If  ie-B'l  «  90°  or  270° ,  go  to  30 

22.  Calculate  m  =  tan  (9-3')/  the  slope  of  the  line  of 
bearing. 

23.  Calculate  c  =  mU+V,  where  (U,V)  is  the  sensor  loca¬ 
tion  for  this  iteration,  c  is  the  y-intercept  of  the  line 
of  bearing. 

2  2  2  2 

24.  Calculate  m  A  +B  -c  .  If  this  quantity  is  less 
than  zero,  the  line  of  bearing  s'  does  not  intersect  the 
threat  ellipse.  Go  to  step  20  and  process  next  increment  of 
size  A3,  if  any  remain.  If  all  I  intervals  of  size  A3  on 
this  side  of  3  have  been  considered  go  to  step  40. 

25.  Calculate  the  points  of  intersection  of  the  line 
of  bearing  with  the  expanded  ellipse: 


(-me  +  (B'/A' ) 


m2  +  B 


'2  -  c2) (A'2) 


A ' 2m2  +  B'2 


(-me  -  (B’/A  ')fh' 

- ;",2  2 

A'  m 


2m2  +  B 
+  B'2 


’2  -  c2)  (A'2) 
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Y1  =  mX^  +  c,  Y2  =  mX2  +  c 

26.  Go  to  32 

27.  If  j V |  >  B * ,  the  bearing  S'  parallel  to  the  x-axis 
does  not  intersect  the  threat  ellipse.  Go  to  step  20  and 
process  next  increment  of  size  AS  if  any  remain.  If  all  -I 
intervals  if  size  AS  on  this  side  of  3  have  been  considered 
go  to  40. 

28.  Calculate  the  points  of  intersection 


29.  Go  to  32 

30.  If  | U !  >  A',  the  bearing  S'  parallel  to  the  y-axis 
does  not  intersect  the  threat  ellipse.  Go  to  step  20  and 
process  the  next  increment  of  size  AS  if  any  remain.  If 
all  I  intervals  of  size  A3  on  this  side  of  3  have  been  con¬ 
sidered,  go  to  40. 

31.  Calculate  the  points  of  intersection: 


32.  Calculate  the  median  width  of  the  strip  of  the 
ellipse  defined  by  the  angular  subinterval  under  considera¬ 
tion.  If  the  subinterval  is  of  size  A3: 

^1  +  ^2  1 

W  =  — ~2 — —(2  tan  ±aB) 

If  the  subinterval  is  of  size  6 A3: 

^1  +  ^2  1 

W  =  — 2 .  (2  tan  ±5A3) 

In  these  expressions  d^  is  the  distance  from  the  sensor  to 
the  ith  point  of  intersection,  i  =  1,  2:  d^  =  ]/{U-X^)2  +  (V-Y^) 

33.  If  j  6-6  * |  =  0°  or  180° ,  go  to  36.  If  j  9-S  * }  =  90° 
or  270°,  go  to  35. 

34.  Evaluate  the  target  density  from  Y. '  to  Y  '  along 

■1  i. 

the  line  y  =  mx  +  c: 

Let  Yl ’  =  min(Y1,Y2) 

Y2'  =  max(Y1,Y2) 

${z)  =  density  function  of  a  standard  normal 

random  variable  evaluated  at  z 

=  (m/foY2+  m2ax2)  $  (c//oy2  +  m2ox2) 

Subdivide  the  interval  (Y^',Y2')  into  n  segments  of  length 
h.  At  the  midpoint  of  each  segement,  y.,  compute: 


when  the  integral  must  be  numerically  evaluated. 
Sum  the  results  over  all  n  segments. 

Go  to  37. 


35.  Evaluate  the  target  density  from  Y^ '  to  ’  along 
the  line  x  =  U: 


Let  Yx'  =  min(Y1,Y2) 


Y2’  =  maxfY-^Y-) 


o(z)  =  density  function  of  standard  normal  random 
variable 


-.  1  ,  U  . 

*■  -  —  $( — } 

°x  °x 


Subdivide  the  interval  (Y^*,^*)  into  n  segments  of  length 


h.  At  the  midpoint  of  each  segment,  y^  ,  compute: 


hK  ~($(-i)) 


'X  -Y 


'j60  .  (str)2cos(e-v)  -  2stTUaos(8-;i<) 

/  - 7—1 - 


X 


2  2 

(stjJ  sin  (s-’i)  -  2stLVjSin(5->')  d  )_ 


where  the  integral  must  be  numerically  evaluated. 
Sum  the  results  over  all  n  segments.  Go  to  37. 


36.  Evaluate  the  target  density  from  X1 '  to  X2  ‘  along 
the  line  y  =  V: 

Let  X^'  =  min(X^,X2) 

X2'  =  maxtX^X  ) 

^(2)  =  density  of  standard  normal  random 
variable 


Subdivide  the  interval  (X^',X2')  into  n  segments  of  length 
h.  At  the  midpoint  of  each  segment,  ,  compute: 


,  x.  360 

hK  ~  (*(-!))  f  exp[- 
“X  CX  6 


1  tStX.^  2qos  ~  2stxxj  (cos 


2  2 

(stjJ  sin  (s-: p)  -  2st£Vsin(S-v)  ^tI. 

+ - 2  •)]  3^ 

"y 


where  the  Integra?  must  be  numerically  evaluated.  Sum 
results  over  all  n  segments. 

37.  Multiply  the  value  of  the  target  density  just 
computed  by  W. 

38.  Multiply  this  result  by 


where  $  is  the  standard  normal  density  function. 


39.  Let  PROS  =  PROB+  (the  results  of  the  calculations 
in  steps  20  through  38  for  each  if  the  I  sufcintervals  of 
size  a S) . 

40.  If  the  fractional  interval  of  size  6 AS  on  this 
side  has  be^n  considered  go  to  41.  Otherwise,  let 

3'  =  jo£3  +  jA3.  Repeat  calculations  21  through  38  once. 
Let  PROS  =  PROB+  (the  result  calculated  at  this  step) . 

41.  If  the  computations  on  both  sides  of  3  have  been 
computed,  go  to  42.  Otherwise  repeat  the  computations 
from  18  to  40  on  the  other  side  of  3  by  letting  as  =  -A3, 

5 AS  =  -oA3. 

42.  Let  3'  =3.  Repeat  steps  21  through  38  once. 

Let  PROB  =  PROB+  (the  results  of  this  calculation) . 

43.  Let  j  =  the  number  of  discrete  threat  speeds  to 
be  considered.  Let  PROB1  =  PROB1  +  P[S=s]PROB.  Go  to  15. 

44.  Repeat  steps  15  through  43  once  with  the  sensor 
located  at  each  of  three  points: 


(1) 

c 

< 

II 

,v  ) ,  the  mean  of  the  sensor 

o 

o 

density; 

(2) 

(U,V)  = 

!uo 

+  ,97ocos (9-3-90) , 

v 

o 

+  .97osin(e-s-90) )  ; 

(3) 

(U,V)  = 

(uo 

-  .97  cos (9-3-90) , 

v  -.97  sin (9-3-90) ) ; 
o 

where  o  is  the  standard  deviation  of  the  sensor  density. 
The  value  of  PROB1  calculated  at  each  iteration  will  be 
weighted  by  the  approximate  probability  that  the  sensor  is 
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located  in  the  region  of  the  sensor  error  circle  repre¬ 
sented  by  the  applicable  value  (U,V) .  If  p^  is  of  the 
order  .86  or  greater,  multiply  by  1/3. 

45.  The  value  of  PR0B2  calculated  after  completion 
of  step  44  is  the  relative  probability  that  the  threat 
lies  along  bearing  3  at  time  t^,  given  the  threat  was  in 
p,  x  100%  ellipse  at  time  t  ,  3.  is  in  (3T,sTJ  and  the 
sensor  is  in  a  p3  *  100%  confidence  circle. 


Fi^ncj3z$ 
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